The potential of historical spy-satellite imagery to support research in ecology and conservation.

Remote sensing data are important for assessing ecological change, but their value is often restricted by their limited temporal coverage. Major historical events that affected the environment, such as those associated with colonial history, World War II, or the Green Revolution are not captured by modern remote sensing. In the present article, we highlight the potential of globally available black-and-white satellite photographs to expand ecological and conservation assessments back to the 1960s and to illuminate ecological concepts such as shifting baselines, time-lag responses, and legacy effects. This historical satellite photography can be used to monitor ecosystem extent and structure, species' populations and habitats, and human pressures on the environment. Even though the data were declassified decades ago, their use in ecology and conservation remains limited. But recent advances in image processing and analysis can now unlock this research resource. We encourage the use of this opportunity to address important ecological and conservation questions.

Legacies of millennial-scale climate oscillations in contemporary biodiversity in eastern North America.

The Atlantic meridional overturning circulation (AMOC) has caused significant climate changes over the past 90 000 years. Prior work has hypothesized that these millennial-scale climate variations effected past and contemporary biodiversity, but the effects are understudied. Moreover, few biogeographic models have accounted for uncertainties in palaeoclimatic simulations of millennial-scale variability. We examine whether refuges from millennial-scale climate oscillations have left detectable legacies in the patterns of contemporary species richness in eastern North America. We analyse 13 palaeoclimate estimates from climate simulations and proxy-based reconstructions as predictors for the contemporary richness of amphibians, passerine birds, mammals, reptiles and trees. Results suggest that past climate changes owing to AMOC variations have left weak but detectable imprints on the contemporary richness of mammals and trees. High temperature stability, precipitation increase, and an apparent climate fulcrum in the southeastern United States across millennial-scale climate oscillations aligns with high biodiversity in the region. These findings support the hypothesis that the southeastern United States may have acted as a biodiversity refuge. However, for some taxa, the strength and direction of palaeoclimate-richness relationships varies among different palaeoclimate estimates, pointing to the importance of palaeoclimatic ensembles and the need for caution when basing biogeographic interpretations on individual palaeoclimate simulations. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

Mapping multiscale breeding bird species distributions across the United States and evaluating their conservation applications.

Species distribution models are vital to management decisions that require understanding habitat use patterns, particularly for species of conservation concern. However, the production of distribution maps for individual species is often hampered by data scarcity, and existing species maps are rarely spatially validated due to limited occurrence data. Furthermore, community-level maps based on stacked species distribution models lack important community assemblage information (e.g., competitive exclusion) relevant to conservation. Thus, multispecies, guild, or community models are often used in conservation practice instead. To address these limitations, we aimed to generate fine-scale, spatially continuous, nationwide maps for species represented in the North American Breeding Bird Survey (BBS) between 1992 and 2019. We developed ensemble models for each species at three spatial resolutions-0.5, 2.5, and 5 km-across the conterminous United States. We also compared species richness patterns from stacked single-species models with those of 19 functional guilds developed using the same data to assess the similarity between predictions. We successfully modeled 192 bird species at 5-km resolution, 160 species at 2.5-km resolution, and 80 species at 0.5-km resolution. However, the species we could model represent only 28%-56% of species found in the conterminous US BBSs across resolutions owing to data limitations. We found that stacked maps and guild maps generally had high correlations across resolutions (median = 84%), but spatial agreement varied regionally by resolution and was most pronounced between the East and West at the 5-km resolution. The spatial differences between our stacked maps and guild maps illustrate the importance of spatial validation in conservation planning. Overall, our species maps are useful for single-species conservation and can support fine-scale decision-making across the United States and support community-level conservation when used in tandem with guild maps. However, there remain data scarcity issues for many species of conservation concern when using the BBS for single-species models.

Rising wildfire risk to houses in the United States, especially in grasslands and shrublands.

Wildfire risks to homes are increasing, especially in the wildland-urban interface (WUI), where wildland vegetation and houses are in close proximity. Notably, we found that more houses are exposed to and destroyed by grassland and shrubland fires than by forest fires in the United States. Destruction was more likely in forest fires, but they burned less WUI. The number of houses within wildfire perimeters has doubled since the 1990s because of both housing growth (47% of additionally exposed houses) and more burned area (53%). Most exposed houses were in the WUI, which grew substantially during the 2010s (2.6 million new WUI houses), albeit not as rapidly as before. Any WUI growth increases wildfire risk to houses though, and more fires increase the risk to existing WUI houses.

Social vulnerability of the people exposed to wildfires in U.S. West Coast states.

Understanding of the vulnerability of populations exposed to wildfires is limited. We used an index from the U.S. Centers for Disease Control and Prevention to assess the social vulnerability of populations exposed to wildfire from 2000-2021 in California, Oregon, and Washington, which accounted for 90% of exposures in the western United States. The number of people exposed to fire from 2000-2010 to 2011-2021 increased substantially, with the largest increase, nearly 250%, for people with high social vulnerability. In Oregon and Washington, a higher percentage of exposed people were highly vulnerable (>40%) than in California (~8%). Increased social vulnerability of populations in burned areas was the primary contributor to increased exposure of the highly vulnerable in California, whereas encroachment of wildfires on vulnerable populations was the primary contributor in Oregon and Washington. Our results emphasize the importance of integrating the vulnerability of at-risk populations in wildfire mitigation and adaptation plans.

The global wildland-urban interface.

The wildland-urban interface (WUI) is where buildings and wildland vegetation meet or intermingle1,2. It is where human-environmental conflicts and risks can be concentrated, including the loss of houses and lives to wildfire, habitat loss and fragmentation and the spread of zoonotic diseases3. However, a global analysis of the WUI has been lacking. Here, we present a global map of the 2020 WUI at 10 m resolution using a globally consistent and validated approach based on remote sensing-derived datasets of building area4 and wildland vegetation5. We show that the WUI is a global phenomenon, identify many previously undocumented WUI hotspots and highlight the wide range of population density, land cover types and biomass levels in different parts of the global WUI. The WUI covers only 4.7% of the land surface but is home to nearly half its population (3.5 billion). The WUI is especially widespread in Europe (15% of the land area) and the temperate broadleaf and mixed forests biome (18%). Of all people living near 2003-2020 wildfires (0.4 billion), two thirds have their home in the WUI, most of them in Africa (150 million). Given that wildfire activity is predicted to increase because of climate change in many regions6, there is a need to understand housing growth and vegetation patterns as drivers of WUI change.

Prioritizing global tall forests toward the 30 × 30 goals.

The Global Deal for Nature sets an ambitious goal to protect 30% of Earth's land and ocean by 2030. The 30 × 30 initiative is a way to allocate conservation resources and extend protection to conserve vulnerable and underprotected ecosystems while reducing carbon emissions to combat climate change. However, most prioritization methods for identifying high-value conservation areas are based on thematic attributes and do not consider vertical habitat structure. Global tall forests represent a rare vertical habitat structure that harbors high species richness in various taxonomic groups and is associated with large amounts of aboveground biomass. Global tall forests should be prioritized when planning global protected areas toward reaching the 30 × 30 goals. We examined the spatial distribution of global tall forests based on the Global Canopy Height 2020 product. We defined global tall forests as areas with the average canopy height above 3 thresholds (20, 25, and 30 m). We quantified the spatial distribution and protection level of global tall forests in high-protection zones, where the 30 × 30 goals are being met or are within reach, and low-protection zones, where there is a low chance of reaching 30 × 30 goals. We quantified the protection level by computing the percentage of global tall forest area protected based on the 2017 World Database on Protected Areas. We also determined the global extent and protection level of undisturbed, mature, tall forests based on the 2020 Global Intact Forest Landscapes mask. In most cases, the percentage of protection decreased as forest height reached the top strata. In the low-protection zones, <30% of forests were protected in almost all tall forest strata. In countries such as Brazil, tall forests had a higher percentage of protection (consistently >30%) compared to forests of lower height, presenting a more effective conservation model than in countries such as the United States, where forest protection was almost uniformly <30% across height strata. Our results show an urgent need to target forest conservation in the greatest height strata, particularly in high-protection areas, where most global tall forests are found. Vegetation vertical structure can inform the decision-making process toward the 30 × 30 goals because it can be used to identify areas of high conservation value for biodiversity protection which also contribute to carbon sequestration.

Priorización de bosques globales altos hacia las metas 30 por 30 Resumen El Tratado Global por la Naturaleza establece una meta ambiciosa de proteger 30% de los continentes y océanos de la Tierra para 2030. La iniciativa 30 por 30 es una forma de asignar recursos para la conservación y extender la protección para conservar ecosistemas vulnerables y sin protección al tiempo que se controlan las emisiones de carbono para combatir el cambio climático. Sin embargo, la mayoría de los métodos de priorización para identificar áreas de elevado valor de conservación se basan en atributos temáticos y no consideran la estructura vertical del hábitat. Los bosques altos globales representan un estructura de hábitat vertical rara que alberga alta riqueza de especies de varios grupos taxonómicos y se asocia con grandes cantidades de biomasa aérea. Los bosques altos globales deberían ser priorizados cuando se planifican áreas protegidas globales en el esfuerzo por alcanzar las metas 30 por 30. Examinamos la distribución espacial de bosques globales con base en el producto Altura de Dosel Global 2020. Definimos a los bosques altos globales como áreas con una altura de dosel promedio por arriba de 3 umbrales (20, 25 y 30 m). Cuantificamos la distribución espacial y el nivel de protección de los bosques altos globales en zonas con gran protección, donde se están alcanzando las metas 30 por 30. Cuantificamos el nivel de protección registrando el porcentaje de bosque alto global protegido con base en la Base de Datos Mundial de Áreas Protegidas 2017. También determinamos la extensión global y el nivel de protección de bosques altos, maduros, no perturbados con base en la mascarilla Paisajes Forestales Globales Intactos 2020. En la mayoría de los casos, el porcentaje de protección decreció a medida que la altura del bosque llegaba al estrato superior. En las zonas poco protegidas, >30% de los bosques estaban protegidos en casi todos los estratos de bosque alto. En países como Brasil, los bosques altos tuvieron un mayor porcentaje de protección (>30% consistentemente) que los bosques de menor altura, presentando un modelo de conservación más efectivo que en países como los Estados Unidos, donde la protección de bosques fue casi uniformemente >30% en los tres estratos de altura. Nuestros resultados muestran una urgente necesidad de enfocar la conservación de bosques en los estratos más altos, particularmente en las áreas muy protegidas, donde se encuentra la mayoría de bosques altos globales. La estructura vertical de la vegetación puede proporcionar información al proceso de toma de decisiones con miras a las metas 30 por 30 debido a que puede ser utilizada para identificar áreas de elevado valor de conservación para la protección de la biodiversidad que también contribuya al secuestro de carbono.

Beyond "greening" and "browning": Trends in grassland ground cover fractions across Eurasia that account for spatial and temporal autocorrelation.

Grassland ecosystems cover up to 40% of the global land area and provide many ecosystem services directly supporting the livelihoods of over 1 billion people. Monitoring long-term changes in grasslands is crucial for food security, biodiversity conservation, achieving Land Degradation Neutrality goals, and modeling the global carbon budget. Although long-term grassland monitoring using remote sensing is extensive, it is typically based on a single vegetation index and does not account for temporal and spatial autocorrelation, which means that some trends are falsely identified while others are missed. Our goal was to analyze trends in grasslands in Eurasia, the largest continuous grassland ecosystems on Earth. To do so, we calculated Cumulative Endmember Fractions (annual sums of monthly ground cover fractions) derived from MODIS 2002-2020 time series, and applied a new statistical approach PARTS that explicitly accounts for temporal and spatial autocorrelation in trends. We examined trends in green vegetation, non-photosynthetic vegetation, and soil ground cover fractions considering their independent change trajectories and relations among fractions over time. We derived temporally uncorrelated pixel-based trend maps and statistically tested whether observed trends could be explained by elevation, land cover, SPEI3, climate, country, and their combinations, all while accounting for spatial autocorrelation. We found no statistical evidence for a decrease in vegetation cover in grasslands in Eurasia. Instead, there was a significant map-level increase in non-photosynthetic vegetation across the region and local increases in green vegetation with a concomitant decrease in soil fraction. Independent environmental variables affected trends significantly, but effects varied by region. Overall, our analyses show in a statistically robust manner that Eurasian grasslands have changed considerably over the past two decades. Our approach enhances remote sensing-based monitoring of trends in grasslands so that underlying processes can be discerned.

Modeled distribution shifts of North American birds over four decades based on suitable climate alone do not predict observed shifts.

As climate change alters the global environment, it is critical to understand the relationship between shifting climate suitability and species distributions. Key questions include whether observed changes in population abundance are aligned with the velocity and direction of shifts predicted by climate suitability models and if the responses are consistent among species with similar ecological traits. We examined the direction and velocity of the observed abundance-based distribution centroids compared with the model-predicted bioclimatic distribution centroids of 250 bird species across the United States from 1969 to 2011. We hypothesized that there is a significant positive correlation in both direction and velocity between the observed and the modeled shifts. We then tested five additional hypotheses that predicted differential shifting velocity based on ecological adaptability and climate change exposure. Contrary to our hypotheses, we found large differences between the observed and modeled shifts among all studied bird species and within specific ecological guilds. However, temperate migrants and habitat generalist species tended to have higher velocity of observed shifts than other species. Neotropical migratory and wetland birds also had significantly different observed velocities than their counterparts, which may be due to their climate change exposure. The velocity based on modeled bioclimatic suitability did not exhibit significant differences among most guilds. Boreal forest birds were the only guild with significantly faster modeled-shifts than the other groups, suggesting an elevated conservation risk for high latitude and altitude species. The highly idiosyncratic species responses to climate and the mismatch between shifts in modeled and observed distribution centroids highlight the challenge of predicting species distribution change based solely on climate suitability and the importance of non-climatic factors traits in shaping species distributions.

Rural land abandonment is too ephemeral to provide major benefits for biodiversity and climate.

Hundreds of millions of hectares of cropland have been abandoned globally since 1950 due to demographic, economic, and environmental changes. This abandonment has been seen as an important opportunity for carbon sequestration and habitat restoration; yet those benefits depend on the persistence of abandonment, which is poorly known. Here, we track abandonment and recultivation at 11 sites across four continents using annual land-cover maps for 1987-2017. We find that abandonment is largely fleeting, lasting on average only 14.22 years (SD = 1.44). At most sites, we project that >50% of abandoned croplands will be recultivated within 30 years, precluding the accumulation of substantial amounts of carbon and biodiversity. Recultivation resulted in 30.84% less abandonment and 35.39% less carbon accumulated by 2017 than expected without recultivation. Unless policymakers take steps to reduce recultivation or provide incentives for regeneration, abandonment will remain a missed opportunity to reduce biodiversity loss and climate change.

Mapping breeding bird species richness at management-relevant resolutions across the United States.

Human activities alter ecosystems everywhere, causing rapid biodiversity loss and biotic homogenization. These losses necessitate coordinated conservation actions guided by biodiversity and species distribution spatial data that cover large areas yet have fine-enough resolution to be management-relevant (i.e., ≤5 km). However, most biodiversity products are too coarse for management or are only available for small areas. Furthermore, many maps generated for biodiversity assessment and conservation do not explicitly quantify the inherent tradeoff between resolution and accuracy when predicting biodiversity patterns. Our goals were to generate predictive models of overall breeding bird species richness and species richness of different guilds based on nine functional or life-history-based traits across the conterminous United States at three resolutions (0.5, 2.5, and 5 km) and quantify the tradeoff between resolution and accuracy and, hence, relevance for management of the resulting biodiversity maps. We summarized 18 years of North American Breeding Bird Survey data (1992-2019) and modeled species richness using random forests, including 66 predictor variables (describing climate, vegetation, geomorphology, and anthropogenic conditions), 20 of which we newly derived. Among the three spatial resolutions, the percentage variance explained ranged from 27% to 60% (median = 54%; mean = 57%) for overall species richness and 12% to 87% (median = 61%; mean = 58%) for our different guilds. Overall species richness and guild-specific species richness were best explained at 5-km resolution using ~24 predictor variables based on percentage variance explained, symmetric mean absolute percentage error, and root mean square error values. However, our 2.5-km-resolution maps were almost as accurate and provided more spatially detailed information, which is why we recommend them for most management applications. Our results represent the first consistent, occurrence-based, and nationwide maps of breeding bird richness with a thorough accuracy assessment that are also spatially detailed enough to inform local management decisions. More broadly, our findings highlight the importance of explicitly considering tradeoffs between resolution and accuracy to create management-relevant biodiversity products for large areas.

Statistical tests for non-independent partitions of large autocorrelated datasets.

Large sets of autocorrelated data are common in fields such as remote sensing and genomics. For example, remote sensing can produce maps of information for millions of pixels, and the information from nearby pixels will likely be spatially autocorrelated. Although there are well-established statistical methods for testing hypotheses using autocorrelated data, these methods become computationally impractical for large datasets. • The method developed here makes it feasible to perform F-tests, likelihood ratio tests, and t-tests for large autocorrelated datasets. The method involves subsetting the dataset into partitions, analyzing each partition separately, and then combining the separate tests to give an overall test. • The separate statistical tests on partitions are non-independent, because the points in different partitions are not independent. Therefore, combining separate analyses of partitions requires accounting for the non-independence of the test statistics among partitions. • The methods can be applied to a wide range of data, including not only purely spatial data but also spatiotemporal data. For spatiotemporal data, it is possible to estimate coefficients from time-series models at different spatial locations and then analyze the spatial distribution of the estimates. The spatial analysis can be simplified by estimating spatial autocorrelation directly from the spatial autocorrelation among time series.

The wildland-urban interface in the United States based on 125 million building locations.

The wildland-urban interface (WUI) is the focus of many important land management issues, such as wildfire, habitat fragmentation, invasive species, and human-wildlife conflicts. Wildfire is an especially critical issue, because housing growth in the WUI increases wildfire ignitions and the number of homes at risk. Identifying the WUI is important for assessing and mitigating impacts of development on wildlands and for protecting homes from natural hazards, but data on housing development for large areas are often coarse. We created new WUI maps for the conterminous United States based on 125 million individual building locations, offering higher spatial precision compared to existing maps based on U.S. census housing data. Building point locations were based on a building footprint data set from Microsoft. We classified WUI across the conterminous United States at 30-m resolution using a circular neighborhood mapping algorithm with a variable radius to determine thresholds of housing density and vegetation cover. We used our maps to (1) determine the total area of the WUI and number of buildings included, (2) assess the sensitivity of WUI area included and spatial pattern of WUI maps to choice of neighborhood size, (3) assess regional differences between building-based WUI maps and census-based WUI maps, and (4) determine how building location accuracy affected WUI map accuracy. Our building-based WUI maps identified 5.6%-18.8% of the conterminous United States as being in the WUI, with larger neighborhoods increasing WUI area but excluding isolated building clusters. Building-based maps identified more WUI area relative to census-based maps for all but the smallest neighborhoods, particularly in the north-central states, and large differences were attributable to high numbers of non-housing structures in rural areas. Overall WUI classification accuracy was 98.0%. For wildfire risk mapping and for general purposes, WUI maps based on the 500-m neighborhood represent the original Federal Register definition of the WUI; these maps include clusters of buildings in and adjacent to wildlands and exclude remote, isolated buildings. Our approach for mapping the WUI offers flexibility and high spatial detail and can be widely applied to take advantage of the growing availability of high-resolution building footprint data sets and classification methods.

Forest phenoclusters for Argentina based on vegetation phenology and climate.

Forest biodiversity conservation and species distribution modeling greatly benefit from broad-scale forest maps depicting tree species or forest types rather than just presence and absence of forest, or coarse classifications. Ideally, such maps would stem from satellite image classification based on abundant field data for both model training and accuracy assessments, but such field data do not exist in many parts of the globe. However, different forest types and tree species differ in their vegetation phenology, offering an opportunity to map and characterize forests based on the seasonal dynamic of vegetation indices and auxiliary data. Our goal was to map and characterize forests based on both land surface phenology and climate patterns, defined here as forest phenoclusters. We applied our methodology in Argentina (2.8 million km2 ), which has a wide variety of forests, from rainforests to cold-temperate forests. We calculated phenology measures after fitting a harmonic curve of the enhanced vegetation index (EVI) time series derived from 30-m Sentinel 2 and Landsat 8 data from 2018-2019. For climate, we calculated land surface temperature (LST) from Band 10 of the thermal infrared sensor (TIRS) of Landsat 8, and precipitation from Worldclim (BIO12). We performed stratified X-means cluster classifications followed by hierarchical clustering. The resulting clusters separated well into 54 forest phenoclusters with unique combinations of vegetation phenology and climate characteristics. The EVI 90th percentile was more important than our climate and other phenology measures in providing separability among different forest phenoclusters. Our results highlight the potential of combining remotely sensed phenology measures and climate data to improve broad-scale forest mapping for different management and conservation goals, capturing functional rather than structural or compositional characteristics between and within tree species. Our approach results in classifications that go beyond simple forest-nonforest in areas where the lack of detailed ecological field data precludes tree species-level classifications, yet conservation needs are high. Our map of forest phenoclusters is a valuable tool for the assessment of natural resources, and the management of the environment at scales relevant for conservation actions.

Early warning sign of forest loss in protected areas.

As humanity is facing the double challenge of species extinctions and climate change, designating parts of forests as protected areas is a key conservation strategy.1-4 Protected areas, encompassing 14.9% of the Earth's land surface and 19% of global forests, can prevent forest loss but do not do so perfectly everywhere.5-12 The reasons why protection only works in some areas are difficult to generalize: older and newer parks, protected areas with higher and lower suitability for agriculture, and more and less strict protection can be more effective at preventing forest loss than their counterparts.6,8,9,12-16 Yet predicting future forest loss within protected areas is crucial to proactive conservation. Here, we identify an early warning sign of subsequent forest loss, based on forest loss patterns in strict protected areas and their surrounding landscape worldwide, from 2000 to 2018.17,18 We found that a low level in the absolute forest cover immediately outside of a protected area signals a high risk of future forest loss inside the protected area itself. When the amount of forest left outside drops to <20%, the protected area is likely to experience rates of forest loss matching those in the wider landscape, regardless of its protection status (e.g., 5% loss outside will be matched by 5% loss inside). This knowledge could be used to direct funding to protected areas threatened by imminent forest loss, helping to proactively bolster protection to prevent forest loss, especially in countries where detailed information is lacking.

Recent collapse of crop belts and declining diversity of US agriculture since 1840.

Over the last century, US agriculture greatly intensified and became industrialized, increasing in inputs and yields while decreasing in total cropland area. In the industrial sector, spatial agglomeration effects are typical, but such changes in the patterns of crop types and diversity would have major implications for the resilience of food systems to global change. Here, we investigate the extent to which agricultural industrialization in the United States was accompanied by agglomeration of crop types, not just overall cropland area, as well as declines in crop diversity. Based on county-level analyses of individual crop land cover area in the conterminous United States from 1840 to 2017, we found a strong and abrupt spatial concentration of most crop types in very recent years. For 13 of the 18 major crops, the widespread belts that characterized early 20th century US agriculture have collapsed, with spatial concentration increasing 15-fold after 2002. The number of counties producing each crop declined from 1940 to 2017 by up to 97%, and their total area declined by up to 98%, despite increasing total production. Concomitantly, the diversity of crop types within counties plummeted: in 1940, 88% of counties grew >10 crops, but only 2% did so in 2017, and combinations of crop types that once characterized entire agricultural regions are lost. Importantly, declining crop diversity with increasing cropland area is a recent phenomenon, suggesting that corresponding environmental effects in agriculturally dominated counties have fundamentally changed. For example, the spatial concentration of agriculture has important consequences for the spread of crop pests, agrochemical use, and climate change. Ultimately, the recent collapse of most agricultural belts and the loss of crop diversity suggest greater vulnerability of US food systems to environmental and economic change, but the spatial concentration of agriculture may also offer environmental benefits in areas that are no longer farmed.

Responses to land cover and grassland management vary across life-history stages for a grassland specialist.

Grassland birds have exhibited dramatic and widespread declines since the mid-20th century. Greater prairie chickens (Tympanuchus cupido pinnatus) are considered an umbrella species for grassland conservation and are frequent targets of management, but their responses to land use and management can be quite variable. We used data collected during 2007-2009 and 2014-2015 to investigate effects of land use and grassland management practices on habitat selection and survival rates of greater prairie chickens in central Wisconsin, USA. We examined habitat, nest-site, and brood-rearing site selection by hens and modeled effects of land cover and management on survival rates of hens, nests, and broods. Prairie chickens consistently selected grassland over other cover types, but selection or avoidance of management practices varied among life-history stages. Hen, nest, and brood survival rates were influenced by different land cover types and management practices. At the landscape scale, hens selected areas where brush and trees had been removed during the previous year, which increased hen survival. Hens selected nest sites in hay fields and brood-rearing sites in burned areas, but prescribed fire had a negative influence on hen survival. Brood survival rates were positively associated with grazing and were highest when home ranges contained ≈15%-20% shrub/tree cover. The effects of landscape composition on nest survival were ambiguous. Collectively, our results highlight the importance of evaluating responses to management efforts across a range of life-history stages and suggest that a variety of management practices are likely necessary to provide structurally heterogeneous, high-quality habitat for greater prairie chickens. Brush and tree removal, grazing, hay cultivation, and prescribed fire may be especially beneficial for prairie chickens in central Wisconsin, but trade-offs among life-history stages and the timing of management practices must be considered carefully.

Habitat heterogeneity captured by 30-m resolution satellite image texture predicts bird richness across the United States.

Species loss is occurring globally at unprecedented rates, and effective conservation planning requires an understanding of landscape characteristics that determine biodiversity patterns. Habitat heterogeneity is an important determinant of species diversity, but is difficult to measure across large areas using field-based methods that are costly and logistically challenging. Satellite image texture analysis offers a cost-effective alternative for quantifying habitat heterogeneity across broad spatial scales. We tested the ability of texture measures derived from 30-m resolution Enhanced Vegetation Index (EVI) data to capture habitat heterogeneity and predict bird species richness across the conterminous United States. We used Landsat 8 satellite imagery from 2013-2017 to derive a suite of texture measures characterizing vegetation heterogeneity. Individual texture measures explained up to 21% of the variance in bird richness patterns in North American Breeding Bird Survey (BBS) data during the same time period. Texture measures were positively related to total breeding bird richness, but this relationship varied among forest, grassland, and shrubland habitat specialists. Multiple texture measures combined with mean EVI explained up to 41% of the variance in total bird richness, and models including EVI-based texture measures explained up to 10% more variance than those that included only EVI. Models that also incorporated topographic and land cover metrics further improved predictive performance, explaining up to 51% of the variance in total bird richness. A texture measure contributed predictive power and characterized landscape features that EVI and forest cover alone could not, even though the latter two were overall more important variables. Our results highlight the potential of texture measures for mapping habitat heterogeneity and species richness patterns across broad spatial extents, especially when used in conjunction with vegetation indices or land cover data. By generating 30-m resolution texture maps and modeling bird richness at a near-continental scale, we expand on previous applications of image texture measures for modeling biodiversity that were either limited in spatial extent or based on coarse-resolution imagery. Incorporating texture measures into broad-scale biodiversity models may advance our understanding of mechanisms underlying species richness patterns and improve predictions of species responses to rapid global change.

Self-perpetuating ecological-evolutionary dynamics in an agricultural host-parasite system.

Ecological and evolutionary processes may become intertwined when they operate on similar time scales. Here we show ecological-evolutionary dynamics between parasitoids and aphids containing heritable symbionts that confer resistance against parasitism. In a large-scale field experiment, we manipulated the aphid's host plant to create ecological conditions that either favoured or disfavoured the parasitoid. The result was rapid evolutionary divergence of aphid resistance between treatment populations. Consistent with ecological-evolutionary dynamics, the resistant aphid populations then had reduced parasitism and increased population growth rates. We fit a model to quantify costs (reduced intrinsic rates of increase) and benefits of resistance. We also performed genetic assays on 5 years of field samples that showed persistent but highly variable frequencies of aphid clones containing protective symbionts; these patterns were consistent with simulations from the model. Our results show (1) rapid evolution that is intertwined with ecological dynamics and (2) variation in selection that prevents traits from becoming fixed, which together generate self-perpetuating ecological-evolutionary dynamics.

Pine plantations and five decades of land use change in central Chile.

The expansion of forest plantations is cause for concern because of their environmental effects, and the loss of native forests and agricultural land. Our goal was to quantify the increase in pine plantation, and concomitant loss of native forests, in central Chile since ca. 1960, and to identify in which settings native forests were lost most rapidly. We analyzed aerial photographs from 1955 and 1961, Landsat images from 1975 and 1998, and Google Earth high-resolution satellite images from 2014. To ensure high classification accuracy, we visually interpreted images for a systematic 3-km grid and assigned each point as either 'pine plantation', 'native forest', 'agricultural-livestock lands', or 'other'. We also calculated latitude, longitude, slope, Euclidean distance to the nearest road and to the nearest pulp mill, and the frequency of land use surrounding each point as potential variables to explain observed land use changes. Pine plantations expansion started even before 1960, when 12% of all points were already pine plantations, was particularly rapid from 1975 (18% of sample points) to 1998 (38%), and stabilized thereafter (37% by 2014). From 1975 to 1998 alone, 40% of native forests were replaced by pine plantations, and agricultural-livestock lands declined by 0.7%, 0.9%, 1% per year before 1975, from 1975 to 1998, and after 1998 respectively. Native forests that were surrounded by pine plantations, were most likely to be converted to plantations, and from 1960 to 1975, also native forests near pulp mills. The probability of change from agricultural-livestock lands to pine plantations was mainly influenced by slope, with most agricultural-livestock lands remaining in areas with low slopes.